Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
鈴土 知明
信学技報, 103(740), p.55 - 60, 2004/03
パーティショニング法を用いて質量保存則が成り立つようにした2次元非同期セルオートマトン(CA)を提案する。このCAは分子動力学系に近い性質を備えており、かつ遺伝的アルゴリズムを用いたルール探索に適している。物理・化学系の創発的空間パターン形成現象のメカニズムの解明やそのシミュレーションに関して、このCAの有用性を議論する。
越智 義浩; Golovkin, I.*; Mancini, R.*; Uschmann, I.*; 砂原 淳*; 西村 博明*; 藤田 和久*; Louis, S.*; 中井 光男*; 白神 宏之*; et al.
Review of Scientific Instruments, 74(3), p.1683 - 1687, 2003/03
被引用回数:10 パーセンタイル:49.35(Instruments & Instrumentation)高強度レーザー光を重水素を充填した燃料球に照射することにより生成されるレーザー爆縮コアプラズマ中の、電子温度,電子密度勾配の時間変化を、時間・空間分解X線分光法により明らかにした。実験においてX線分光ストリークラスターカメラによる時間分解X線スペクトルデータと単色X線駒取りカメラによる特定のラインX線に対応する単色のX線エネルギーの時間分解二次元画像データを同時に取得した。これらのデータをもとに遺伝的アルゴリズムを用いた自己無頓着な解析手法により、爆縮コア内部の電子温度,電子密度勾配を求めることに成功した。得られた成果と流体コードによるシミュレーションとを比較し、爆縮コアプラズマ形成の動的過程の診断を進めている。
渡辺 俊樹*; 真田 佳典*; 藪内 聡
JNC TY7400 2000-001, 72 Pages, 2000/03
本研究では、地下に電磁波エネルギーを送り込み、それに対する地下の応答を観測することにより地下構造を解析する地球物理学的調査において、正弦波的に変化する電磁エネルギーの周波数を高精度に制御し、連続して長時間地下に送り込む手法について研究開発を行った。本手法を採用することによって、調査深度の拡大や解析精度の向上ならびに調査機器の小型化を目指した。本報告書は本研究テーマのうち数値シミュレーションおよびデータ解析技術に関する研究の報告である。電磁波動現象は周波数によって挙動が大きく異なるため、周波数が数HZ数KHZ程度の低周波数を利用した電磁探査法および周波数が数十MHZ数GHZ程度の高周波数を利用した地下レーダー法のそれぞれについて電磁波動現象のシミュレーションを行った。次に、本調査手法の主な適用対象を岩盤内のフラクチャおよび含水破砕帯と想定し、地下レーダー法による岩盤内の微細構造の検出の可能性について検討した。次に、計測した電磁波動データから岩盤内部の物性値分布を求めるインバージョン解析手法について、数値実験および実データへの適用を行った。さらに、高分解能な時間-周波数解析アルゴリズムについて検討した。
山田 澄*
JNC TJ8400 2000-051, 122 Pages, 2000/02
使用済み核燃料再処理工場等における安全性の確保のために、安価で即応性の良い臨界安全監視システムの開発が望まれており、基礎研究による知見を基に具体的な臨界安全監視システムとして実現することを目的として開発研究が進められている。本研究に用いている基本的手法は、未臨界体系で観測される中性子信号揺らぎをオンライン計算処理して未臨界度を推定しようとする試みであり、オンラインARMAモデル同定手法を安全監視システムに応用のためのシステム構成要素モジュールの開発、核燃料サイクル開発機構のDCA実験データ解析など実用化に必要となる知見の蓄積も行ってきたが、この間、Microsoft Windowsを基本ソフトとするDOS/V系パーソナル計算機の爆発的な普及により、NECを中心とするMS-DOS系の計算機がその地位を失い、今やMS-DOS系の計算機はその姿を消しつつある。本研究ではMS-DOS系計算機を用いたプログラム開発を行ってきており、その研究成果を今後も有効に利用するためには、これらを商業主義に左右されないOSに基礎を置くワークステーションに移すことが重要な課題となってきた。また、本研究で提案している臨界安全監視システムの設計基本思想は高信頼性、耐雑音性、そして高拡張性にある。近年エンジニアリングワークステーションの性能は飛躍的に向上し、且つ価格もかなり安くなっており、この様にめざましい発展を遂げている計算機に着目すると、信号解析、表示、警報信号発生に、それらのアルゴリズムをROM化したICチップを用いるよりも、高性能の計算機を用意し、各種の信号解析ソフトウェアや、必要に応じて他のシステム同定アルゴリズムを組み合わせた総合的信号処理システムとすることが望ましい。この様な観点から、臨界安全監視用プログラムをワークステーションに移植することは非常に重要である。そこで、今年度より、これまでMS-DOS系計算機上で開発された臨界安全監視用プログラムを順次年次計画を立ててワークステーション上に移植することとなった。今年度は、臨界安全監視システムの中心的ソフトウェアである逐次型ARMAモデル同定プログラム(ARMAX1)の移植作業を行った。なお、本報告書として自己完結性を保つため、理論的な未臨界原子炉のARMAモデル導出を示すとともに、逐次型ARMAモデル同定アルゴリズムに関する要点を再
飯島 隆; 白鳥 芳武; 松本 光雄; 川島 仁*
JNC TN3410 2000-002, 93 Pages, 2000/01
ふげん発電所は新型転換炉の原型炉であり、熱中性子炉におけるプルトニウム利用において、昭和54年3月の運開以降、これまでに各種の混合酸化物(MOX:Mixed Oxide)燃料集合体を原子炉に装荷するなど、核燃料の多様性を実証してきている。また、運転開始以来、燃料集合体の漏えいは1体もなく高い信頼性を得ており、MOX燃料集合体も700体を超える使用実績を有している。この数は熱中性子炉として、世界最大の使用体数を誇っている。しかしながら、新型転換炉開発についてはその役割が終了しつつあることから基本的に撤退し、「ふげん」については平成15年に運転を停止することが決定されている。そのため、限られた運転期間において、過去の技術開発成果を含め、プルトニウム利用技術やプラント管理技術についてとりまとめたものである。
陳 振茂; 青砥 紀身; 加藤 章一
JNC TN9400 99-061, 32 Pages, 1999/07
本報告書では、自然磁束漏洩からき裂・損傷を非破壊的に検査する研究の一環として、測定した磁束信号より材料における磁荷分布(損傷による)の再構成を行った。この代表的な非適切問題には最小自乗法に基づいた反復計算アルゴリズムを用いた。問題を適切化するために初期値、重み係数及び反復計算の回数の選び方を検討した。シミュレーション信号を用いた再構成結果より、本手法がノイズの少ない信号に対して有効であることを確認した。ノイズに対するロバースト性を向上するために、ウェーブレットをガラキン法に適用した手法をシステム方程式の離散化に導入した。最小自乗法と比較した結果、ウェーブレットを用いた手法はS/N比の低い信号に対しても有効であることが判った。本報告書では最小自乗法に基づいた手法を1次元及び2次元の磁荷分布、ウェーブレットを用いた手法を1次元の磁荷の再構成に適用し、提案した手法の妥当性を実証した。
守田 幸路; 飛田 吉春; 近藤 悟; E.A.Fischer*
JNC TN9400 2000-005, 57 Pages, 1999/05
高速炉安全解析コードSIMMER-IIIで使用する解析的状態方程式(EOS)モデルを開発した。汎用的な熱力学的関数式を使用した本モデルは、計算効率を犠牲にすることなく、幅広い温度および圧力領域での炉心物質の熱力学的特性を充分な精度で記述し、基本的な熱力学的関係を満足するように設計されている。本報告書では、このEOSモデルと結合した圧力反復計算の流体力学アルゴリズムについても記述した。二酸化ウラン、混合酸化物燃料、ステンレス鋼およびナトリウムの臨界点までのEOSデータについては、最新でかつ最も信頼できるデータに基づき、基本的な熱力学的関係を用いて求めた。EOSデータの熟力学的整合性と精度についても既存データと比較することで議論した。
Dershowitz, W.*; Doe, T.*; Shuttle, D.*; Eiben, T.*; Fox, A.*; Emsley, S.*; Ahlstrom, E.*
JNC TJ1400 99-012, 582 Pages, 1999/02
本報告書は、亀裂流動モデルの開発と性能評価への適用性に関する平成10年度のGolder Associates Inc.の研究結果を報告するものである。平成10年度のGolder Associatesの研究目的は、第2次取りまとめ(性能評価報告書)に理論的及びモデリングについて支援することである。加えて、Golder AssociatesはJNCのAspoプロジェクトについて技術的にサポートする。性能評価の支援のための主な研究は、流れと移行のシミユレーション、移行経路の簡素化の解析、掘削影響領域による影響の調査、ソフトウェアの検証及び比較検証、モンテカルロシミユレーションの信頼性の解析である。また、フィックの拡散のアルゴリズムを、Laplace Transform Galerkin solute transportに適用したアルゴリズム改良を施した。Aspoプロジェクトに対する支援項目は、TRUE-1での岩石ブロック中での収着性トレーサの移行予測のためのモデリング、TASK5での1kmスケールの地化学的な移行経路解析、TRUEブロックスケール試験のデータ解析と試験計画である。技術的情報については、本報告書に示した。
not registered
JNC TJ1400 99-004, 110 Pages, 1999/02
臨界安全監視システムは高い信頼性と新しい手法を容易に取り込むことのできる柔軟性が必要である。この様な観点から1997年にはシステムの基本構成を報告した。ここでは、この基本デザインに新たに加えるべきシステムと手法について報告すると共に、1996年のDCA実験データの解析結果を報告した。1:本報告書では、先ず、ガンマ線信号を用いたFeynman-法による未臨界度推定モジュールと、その信号取得に必要となる高速パルス信号処理系の追加を提案した。ガンマ線信号に対してFeynma-法を用いることにより原子炉の減衰定数を推定することが可能であることが最近報告されており、このモジュールを追加することは、計数率の広い範囲に亙って臨界安全監視システムの信頼性を高めるものと考えられる。2:昨年度、ARMAモデル同定と平行して信号解析を行う冗長系として最小2乗法に基づくアダプティブフィルタ(ADF)アルゴリズムモジュールを報告したが、この手法は有色雑音が混入する信号に対してパラメータの収束が遅くなることが知られており、ここでは近年注目されているブロックアルゴリズムの基礎理論と、これを実現するための跳躍アルゴリズムについて報告した。3:アナログ信号系では不可避の電源雑音を除去するノッチフィルタを設計法を実際のソフトウェアにつ臨界安全監視システムは高い信頼性と新しい手法を容易に取り込むことのできる柔軟性が必要である。この様な観点から1997年にはシステムの基本構成を報告した。ここでは、この基本デザインに新たに加えるべきシステムと手法について報告すると共に、1996年のDCA実験データの解析結果を報告した。1:本報告書では、先ず、ガンマ線信号を用いたFeynman-法による未臨界度推定モジュールと、その信号取得に必要となる高速パルス信号処理系の追加を提案した。ガンマ線信号に対してFeynma-法を用いることにより原子炉の減衰定数を推定することが可能であることが最近報告されており、このモジュールを追加することは、計数率の広い範囲に亙って臨界安全監視システムの信頼性を高めるものと考えられる。2:昨年度、ARMAモデル同定と平行して信号解析を行う冗長系として最小2乗法に基づくアダプティブフィルタ(ADF)アルゴリズムモジュールを報告したが、この手法は有色雑音が混入する信号に対してパラメータの
Yang Jin An*; 三原 隆嗣
JNC TN9400 99-013, 89 Pages, 1998/12
成功基準が時間と共に変化するシステムの信頼度をより現実的に評価することを目的として、使命時間を複数のフェイズに分割して評価を行うフェイズドミッション解析コード:PHAMMONの開発を行っている。大規模なシステムモデルにも適用可能とするためモンテカルロ法を採用しており、既に強制遷移法と故障バイアス法という2種類の分散低減法を取り入れている。しかしながら、評価対象によってはこれらの方法のみでは分散低減の度合が不十分な場合もあり、本研究ではさらなる改良を目的として、頂上事象成立バイアス法を適用した計算アルゴリズムの改良を行った。頂上事象成立バイアス法では、任意のシステム状態から各々のカットセット成立状態までの遷移の起こりやすさを指標化した「遷移距離」を計算し、機器の運転成功から故障状態への状態遷移のサンプリングを「遷移距離」の最も短いカットセット成立に向けてバイアス(偏向)させることにより、カットセットが成立する有効ヒストリーをより効率的に発生させようとするものである。上述の頂上事象成立バイアス法をPHAMMONコードに導入し、大型高速増殖炉モデルプラントの崩壊熱除去系を評価例として適用計算を実施した。その結果、本方法はモンテカルロ法による出力結果の分散をさらに低減する上で有効であるとの結論を得た。
not registered
PNC TJ1632 98-001, 112 Pages, 1998/03
安価で即応性の良い臨界安全監視システムの開発が望まれている。昨年度は、過去6年間の研究によって得られた知見を基に、実際の臨界安全監視システムとしてシステム化するために必要となるシステムの構成要素モジュールに関する検討をおこなった。そこで、本年度は昨年度提案した解析表示モジュールの中で、ARMAモデル同定と平行して信号解析を行う冗長系として重要な要素となるアダプティブフィルタ(ADF)アルゴリズムモジュールの開発と、このモジュールを用いた実データ解析を行い、ADFアルゴリズムの未臨界度推定手法としての有効性を検討し、以下の成果を得た。ADFによるシステム同定においては、修正係数が大きい時には推定されるパラメータは入力データの統計的性質の変化に素早く追従するが、その反面推定されたパラメータの統計的変動が大きい。逆に修正係数を小さくすると統計的変動は小さくなるが、実際にシステムパラメータに変動がある場合には追従が遅くなり、この現象は定常時系列、非定常時系列を問わず、全てのADFパラメータ推定結果について言える。推定されたパラメータから得られる未臨界度はいずれの場合においても統計的変動が非常に大きい。この統計的変動は逐次型ARMAモデル同定アルゴリズムによる未臨界度推定の場合に比べて非常に大きく、ADFアルゴリズムによる結果をそのまま未臨界度推定値とするにはやや問題がある。従って、推定結果を平滑化する二次的なアルゴリズムを付加することが重要になる。本研究では、ADFによる推定パラメータに500次の単純移動平均、更に、その平滑化されたパラメータの再移動平均、500データ毎の相加平均などの平滑化を行った結果、いずれも統計的変動を充分小さくすることができた。
椹木 哲夫*
PNC TJ9604 98-001, 83 Pages, 1997/03
現在までに原子力プラントで発生したトラブルの多くは、運転操作ミスなどの人的要因に関係したものである。そこで安全性、信頼性向上を目的として、従来人間の運転員や保守員の果たしてきた役割を人工知能技術を主体とした自動化システムに代替し、できるかぎり人的要因を排除していこうとする自律型プラントの開発研究が盛んに進められている。自律型プラントにおいては、プラントに関する膨大な量の知識を実施し、この知識を実時間で利用可能にしていくものであるが、そこで問題になるのが、これら膨大な知識の獲得と管理である。すなわち合理的な知識として実装できる知識の量には自ずと限界があり、このような不完全な知識を用いる限り、想定外事象の生起に対しては極めて脆弱なものにならざるを得ない。もちろん最終的には人間オペレータの判断に委ねることが考えられるが近年の人間オペレータへの監視業務一辺倒に偏りがちな自動化システムとの役割分担やモニターすべき情報の飛躍的な拡大は自律機能のみならず、オペレータに対して想定外事象が生じた差異の緊急時に冷静沈着な判断を支援することも併ねた知識化が必須となる。すなわち自律型プラントに求められるのは、このような状況で生起している事象を自身の知識の不備をも考慮した解釈を行うことで的確に現時点での矛盾を検知し、その解消策の提案をオペレータに対して提示したり、自らその解消を試みることで、自身の知識を常に更新していける機構の実現である。さらに知識ベースの実時間運用を考える際には、時間的に不変な知識ばかりを想定することは難しく、対象システムの時間的推移に合せて適用すべき知識もダイナミックに追従しながら更新していくことが求められるが、この意味においても、いつどの時点で現状の知識の確認し、これを改定していけばよいかについてのメンタルモデル判断機能が不可欠となる。本報告ではこのような状況に際して、想定外の事象に対して現在の知識ベースの知識との矛盾を速やかに検知し、自律プラント自らが、あるいはオペレータがその矛盾を解消を的確に行うことを可能にするための方策について、人工知能における機械学習(machine learning)分野から、帰納学習、演繹学習を中心に概説する。そしてそれらの学習手法における矛盾検知とその解消による知識獲得がそのように実現されるかについて、1)既存知識の矛盾検知、2)合理的
岡本 孝司*; 鶴 大悟*; 文沢 元雄
Proc. of 1997 ASME Fluids Engineering Division Summer Meeting (FEDSM'97), 0, p.1 - 6, 1997/00
高温ガス炉配管破断時の密度差流れ及び密度の異なる可燃性気体の熱流動を調べる基礎的研究として、3次元密度分布の再構築を行った。本手法は、マッハツェンダー干渉計で作成した2次元干渉縞可視画像を遺伝アルゴリズムを用いた解析プログラムで3次元密度分布を再構築することである。従来の手法である最大エントロピー技術では再構築に12画面の画像データが必要であったが、今回の遺伝アルゴリズムを用いた手法では、3画面の画像データで再構築が可能であることを示した。したがって、本手法は限定された画像データより3次元密度分布を再構築することに有効であることが分かった。
芳野 隆治; Koga, J. K.*; 竹田 辰興
Fusion Technology, 30(2), p.237 - 250, 1996/11
ディスラプションにおけるプラズマ電流消滅時には、真空容器に渦電流が誘起され、プラズマ電流中心の垂直位置(Z)の検出に誤りを発生する。この誤りは、垂直方向のプラズマの移動(VDE)をひきおこし、真空容器内部の損傷の原因となる。このため、ITER等の核融合炉では、VDEの回避が必要とされている。一方、定常状態では、高い精度で、Zを検出しプラズマの高い性能を維持する必要がある。このため、プラズマ電流消滅時と定常状態の両方に適用できるZの検出アルゴリズムについて検討した。この結果、統計処理法を用いると、Zの検出精度を両方の状態にて同時に改善するのが困難であるのに対し、ニューラルネットワーク法は、同時に改善することが可能であることを解明した。
岸本 牧; 坂佐井 馨; 荒 克之; 藤田 隆明; 鈴木 康夫*
IEEE Transactions on Plasma Science, 24(2), p.528 - 538, 1996/04
被引用回数:1 パーセンタイル:4.38(Physics, Fluids & Plasmas)トカマクにおける磁界計測は、内部インダクタンス等を求めるために必要なプラズマ電流分布に関する重要な情報を与える。ところが磁界計測のみからプラズマ電流分布再構成をするのは、解が一意に決まらない問題である。本研究では、磁界計測データからのトカマクプラズマ電流分布構成問題を、真空容器断面への電流最適配置化問題とみなし、この組み合せ最適化問題を解く手法として、遺伝的アルゴリズムを用いた手法とホップフィールド型ニューラルネットワークと遺伝的アルゴリズムの結合手法の2つの方法を用いた。そしてこれらの手法をトカマク型核融合装置JT-60Uのダイバータプラズマに適用し、従来プラズマ電流プロファイル再構成手法として一般的に用いられている平衡コードによって得られたものと非常に良い一致が得られた。
鈴木 篤之*; 長崎 晋也*
PNC TJ1602 96-001, 67 Pages, 1996/03
天然バリアを構成する岩盤の中で人工バリアに交差した亀裂では、亀裂中のある領域まではベントナイトで充填されているものと考えられている。人工バリア中での放射性核種の拡散は亀裂充填ベントナイト中を移行し、亀裂ネットワーク中の地下水流に放出されるが、こうした核種移行プロセスは人工バリアと天然バリアとの間の掘削影響領域に影響されることが予想される。しかし、このプロセスは非常に複雑であるため、掘削影響領域を模擬しアクチニド元素の吸着やコロイド形成機構等を考慮した実験並びにモデルに基づいた数値解析が必要である。そこで本研究では、ベントナイトで充填された仮想的な岩盤亀裂中での種々の酸化還元電位(Eh)に対して、アメリシウム、ネプツニウムの吸着分配係数を実験により測定した。このとき、Ehの関数として実験的に求められたKdを用いることにより、ベントナイトで満たされた亀裂中での放射性核種の移行を数値的に評価することができる。ベントナイトで満たされた亀裂中から天然バリア内亀裂ネットワークへの物質輸送についてシミュレーション解析を行った結果、ウランの同位体と娘核種の移行には影響がないことがわかった。さらに、遺伝的アルゴリズムとファジー理論との融合による解析手法を用いて人工バリアにおける放射性核種の移行に関する不確実性解析を行った。
岸本 牧; 坂佐井 馨; 荒 克之
Journal of Applied Physics, 79(1), p.1 - 7, 1996/01
被引用回数:10 パーセンタイル:48.61(Physics, Applied)磁界逆問題は、本質的に解が一意に定まらないill-posedな問題であり、解析的な手法では解くことが出来ない。そこでこの磁界逆問題を、磁界源電流の存在する空間領域への電流最適配置化問題へと変換し、その最適化問題を解く手法として、ホップフィールド型ニューラルネットワークと遺伝的アルゴリズムを組み合わせた新しい手法を提案した。そしてコンピュータによるシミュレーションによりその有効性を示した。
Ugolini; 吉川 信治; 小澤 健二
PNC TN9410 95-210, 11 Pages, 1995/09
蒸気発生器の出口蒸気温度の正確な制御は、原子力発電所の水/蒸気系全体の性能を向上させるために重要である。本報告書では、制御対象の挙動予測及び制御対象の上流で時間的に先行して観測される境界条件をニューラルネットワークによって制御信号に反映させるアルゴリズムについて述べる。このアルゴリズムはニューラルネットワークを用いたモデル適応制御法(MRACnn)に組み込まれている。MRACnnが原子炉の蒸気発生器のような非線形な機器をも制御できることは既に報じられているが、MRACnnシステムの性能向上のためにニューラルネットワークの特性を更に利用する手法については今まで考えられてはいなかった。本報告書に述べるMRACnnの改良型は高速炉プラントシュミレータと接続して、蒸気発生器出口蒸気温度の制御に用いられた。この結果、制御対象の挙動予測及び制御対象の上流で時間的に先行して観測される境界条件を制御信号に反映させることが制御性能を向上させることを確認した。
not registered
PNC TJ1632 95-001, 107 Pages, 1995/03
安価で即応性の良い臨界安全監視システムの開発が望まれている。本研究は未臨界体系で観測される中性子信号揺らぎを用いてオンラインで未臨界度を推定しようとする試みであり、そのための基礎的研究を平成3年度より行ってきた。未臨界度が浅い状態から深くなる非定常状態では、この逆の非定常状態のときほど高速には未臨界度の変化を追跡できない問題があり、本年度はこの問題を解決するための最適アルゴリズムの検討を行うと共に、動燃事業団DCA実験炉雑音データを解析し、以下の結果を得た。1.未臨界度が変動する非定常状態における未臨界度推定アルゴリズムとしては、我々が提案してきたP行列活性化(P-matrix Activation:PMA)法を更に改良したP行列対角要素活性化法(PMDA法:P-Matrix Diagonal Activation Method)が有効である。特に、未臨界度が深い方向に変化する場合は最小2乗関数を評価関数とするELS法、及び適切な忘却計数パラメータを用いたRPE法が優れていると言える。しかし、未臨界度の変化がどの方向に変化するのか分からないような一般的な場合、細かいパラメータを選択する必要のない点で、RPE法よりもELS法にP行列対角要素活性化を組み合わせたアルゴリズムが最も有効であるととの結論を得た。2.従来から我々が行ってきた近大炉実験データの信号処理法が、DCA実験装置を用いた中性子検出器揺らぎ信号の解析にも有効であることが実証された。すなわち、本年度、動燃事業団・大洗工学センターのDCA実験装置を用いて行われた未臨界定常運転時の中性子信号揺らぎデータを解析した結果、未臨界度が0$11.2に対し、周波数スペクトル折れ点周波数は数Hzから数十Hzの間に分布している。
根山 敦史*; 中川 浩一*; 明石 浩一*; 白川 敏彦*
PNC TJ1458 95-004, 22 Pages, 1995/03
本研究は、TRU廃棄物の処分システムの挙動を長期にわたって予測するための個別現象のうち、化学反応と物質移行現象を結合した評価コードを対象として、以下の研究を実施した。(1)地球化学モジュールの機能において、溶解・沈澱が取り扱える様に改良すると共に、1種類の鉱物生成/2種類の平衡計算が可能となる様コードの改良を行った。(2)改良したコードのアルゴリズムの妥当性ならびに数値的な精度を確認するため、検証解析を行った。物質移行、地球化学、これらを結合したモジュールの検証解析は、HYDROGEOCHEM,PHREEQEとのベンチマーク解析を行った。(3)改良したコードの実現象への適用性を確認するため、確証解析を行った。確証解析は、albite/純水、コンクリート/純水のカラム試験データとの比較を行い、改良したコードの適用性について検討した。(4)確証解析による検討結果に基づき、さらに改良すべき点を明確にした。